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ABSTRACT 
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1 r*| Context. This work was conducted as part of the SPECPDR program, dedicated to the study of very small particles and astrochemistry, in 
Cln' Photo-Dissociation Regions (PDRs). 

Aims. We present the analysis of the mid-IR spectro-imagery observations of Ced 201, NCG 7023 East and North-West and p Ophiuchi West 
filament. 

^ Methods. Using the data from all four modules of the InfraRed Spectrograph onboard the Spitzer Space Telescope, we produced a spectral cube 
, ranging from 5 to 35 pm, for each one of the observed PDRs. The resulting cubes were analysed using Blind Signal Separation methods (NMF 

■ and FastICA). 

. Results. For Ced 201, p Ophiuchi West filament and NGC 7023 East, we find that two signals can be extracted from the original data cubes, 

■ which are 5 to 35 pm spectra. The main features of the first spectrum are a strong continuum emission at long wavelengths, and a broad 7.8 
5_l ' pm band. On the contrary, the second spectrum exhibits the classical Aromatic Infrared Bands (AIBs) and no continuum. The reconstructed 

spatial distribution maps show that the latter spectrum is mainly present at the cloud surface, close to the star whereas the first one is located 
slightly deeper inside the PDR. The study of the spectral energy distribution of Ced 201 up to 100 pm suggests that, in cool PDRs, the 5-25 pm 
continuum is carried by Very Small Grains (VSGs). The AIB spectra in the observed objects can be interpreted as the contribution of neutral 
and positively-charged Polycyclic Aromatic Hydrocarbons (PAHs). 

Conclusions. We extracted the 5 to 25 pm emission spectrum of VSGs in cool PDRs, these grains being most likely carbonaceous. We show 
that the variations of the mid-IR (5-35 pm) spectra of PDRs can be explained by the photo-chemical processing of these VSGs and PAHs, 
VSGs being the progenitors of free PAHs. 

Key words, astrochemistry — ISM: dust — ISM: lines and bands — reflection nebulae — infrared: ISM — methods: numerical — methods: 
observational 



o 



1. Introduction 

The Aromatic Infrared Bands (AIBs), observed in emission 
in UV irradiated interstellar matter, have been attributed to 



Polycyclic Aromat ic Hydroca r bons (PAHs) bv lLeger & P uget 



rolycy 

dl984 and lAllamandola et all dl985l) . La ter observations with 
the Infrared Space Observatory (ISO) dKessler et all []996) 
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have largely contributed to the study of these macromolecules 
and their evolution in the ISM, together with many fundamen- 
tal studies involving both laboratory experiments and quan- 
tum chemistry calculations. Although these studies seem to 
confirm that the AIB carriers are PAH-like species, a precise 
identification of the mixture of the interstellar molecules has 
not been possible y e t. In p arallel to this quest for identifica- 



tion, ICesarskv et alj d2000l) gave strong evidence of the exis- 
tence of Very Small carbonaceous Grains (VSGs), from the 
ISO observations of the reflection nebula Ced 201. They pro- 
posed that these grains, emitting broad bands and a continuum, 
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are destroyed by UV radiation and/or shock waves freeing the 
AIB carriers in the process. In addition, they suggested that 
these grains could be present everywhere in the ISM, but de- 
tectable in the mid-IR only under special conditions. Photo- 
Dissociation Regions (PDRs) constitute the transition regions 
between UV-rich regions and molecular clouds. In PDRs with 
limited excitation, the infrared emission is dominated by dust 
particles. To emit at mid-IR wavelengths, these particles must 
be heated to high temperatures implying that they are very 
small (from ~ 1 to 10 nm) and not in thermal equilibrium with 
the incident radiation field. PDRs are extremely useful to probe 
the properties of very small dust particles, because the local 
UV radiation field varies as a function of the distance from 
the star. This change in the excitation conditions is expected 
to imply changes in the dust processing and therefore pro- 
voke significant variations in the observed infrared spectrum. 
The analysis of the observed spectral variations can provide 
valuable clues on the physical and chemical properties of the 
emitting very small particles. The study of such regions with 
a high spectral and spatial resolution was first made possible 
with ISOCAM onboa rd the Infrared Space Observatory (ISO) 
dCesarsky et al.lll996h . Most analyses of the ISOCAM spectral 
cubes have been achieved with simple mathematical decom- 
positi ons of the observed spectra in bands (the AIBs) and con- 
tinua dAbergel et al .120031) . However, more sophisticated meth- 
ods were used elsewhere to describe the observed spectrum on 
each pixel as the combination of the emission of several chem- 
ical populations. This approach requires the use of methods 
such as Single Value Decomposition (SVD is applied to the 
covariance matrix of the data, which is equivalent to Principal 
Component An alysis). Apply i ng SV D to the observations of 



differe nt PDRs, Boissel et al. ( 2001 ) and then iRapacioli et al 



20051) (RJB hereafter) were able to extract the spectra of dif- 
ferent types of very small dust particles: PAH-like species and 
VSGs. Although pioneering, this work was restricted to two 
objects of the ISOCAM database, because of the limited effi- 
ciency of the method. 

In the field of signal processing, quite a large number of 
algorithms have been developed to extract the original signals 
from observations which are mixtures of these signals. Such 
algorithms are referred to as Blind Signal Separation (BSS) 
methods. The most widely used are Independent Component 
Analysis (ICA) approaches, with algorithms like FastICA 
(IHyvarinenl 1 19991) . The use of ICA is so widespread that 
it is sometimes presented as the only class of BSS meth- 
ods, although some oth er algorithms like N on-negative Matrix 
Factorisation (NMF) dLee & Seung 1 19991) have proven to be 
efficient. 

This paper presents the results of a study conducted with 
the newly available spectro-imagery data from the Infrared 
Spectrograph (IRS) (IHouck et al.l 120041) onboard Spitzer. Our 
goal is to use both the spectral and spatial information avail- 
able in IRS data cubes to reveal the origin of the observed 
spectral variations. The IRS spectrograph provides several im- 
provements when compared to ISOCAM. These are a higher 
sensitivity (about 0.06 mJy from 6 to 15 pm, and 0.4 mJy from 
15 to 38 pm), a better spectral resolution (-^j ~ 60 - 127 in 
the low resolution mode) and a range of wavelengths extended 



to 35 pm. We use two classes of BSS algorithms, i.e. ICA and 
NMF, to extract the spectra of the different populations of very 
small particles. This work, carried out on four PDRs, extends 
previous studies on the following aspects: 

i) we have gained in sensibility and spectral resolution; 

ii) a wider range of wavelengths is available; 

iii) the BSS methods are more powerful than S VD-based meth- 
ods; 

iv) as a consequence, our results involve fainter objects. 

After describing the observations and data reduction in 
Sect. 2, we define the BSS methods we have used for our anal- 
ysis in Sect. 3. We then present the results we have obtained 
applying these methods to the data (Sects. 4-5) and provide ev- 
idence of their efficiency (Sect. 6). We show that the extracted 
spectra can be attributed to different populations of dust (Sect. 
7). Using the extracted in the studied PDRs in Sect. 8. 

2. Observations 



2.1. Selected PDRs 

This study was co nducted as part of the SPECPD$\ program 
( Joblin et alj|2005l) in which a total of 1 1 PDRs were observed 
with all the instruments onboard Spitzer. The PDRs selected 
for this paper, Ced 201, the p-Ophiuchi filament, NGC 7023 
East a nd North, hav e mild radiation fields (200 to 1300 in units 
of the iHabingld 19681) . G =1.6 x lO" 10 W cm 2 ) and relatively 
simple geometries. They were chosen among the 1 1 observa- 
tions because their IR emissions show the strongest spectral 
changes across the PDR, thus providing the best information 
on the nature of the carriers. 



2.1.1. Ced 201 

Ced 201 is a compact reflection nebula lying in the Cepheus 
constellation about 420 pc from the Sun (lCasevlll99ll) . Unlike 
common reflection nebulae, which give birth to their own star, 
Ced 201 is the result of a chance encounter between a B9.5 V 
star (BD+69°1231) and a molecular c loud, at arelative velocity 
of about 12 km s 1 rtwitt et al.lll987l) . Because of this particu- 
larity, Ced 201 is a very a ttractive object f or the spectroscopic 
study of interstellar dust. ICesarsky et al. I d2000l) revealed the 
presence of very small carbonaceous grains, with a spectrum 
exhibiting bands and continuum emission, identified as the gen- 
eral mid-IR VSG spectrum by RJB. 

2.1 .2. NGC 7023 East and North-West 

NGC 7023 is a well studied reflection neb ula located at 430 pc 
from the sun dvan den Ancker et al.ll997h . Three PDRs are vis- 
ible around the Herbig Be illuminating star, HD 200775. The 
brightest one, NGC 7023 North-West (NGC 7023-NW here- 
after), lies about 40 " North-West of the star and another one 
lies around 70 " South. The last one, NGC 7023 East (NGC 
7023-E hereafter), is located further East of the star (~170")- 



http://www.cesr.fr/joblin/SPECPDR_public/Home 
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The spectroscopic studies of the North-W est and South PDRs 
we re performed by Cesarskv et al. ( 19961) using ISOCAM, and 
by IWerner et al.l (120041) with IRS, leading to the discovery of 
new infrared features. In the present paper, we analyse the IRS 
spectro-imagery observations of NGC 7023-E and NGC 7023- 
NW. 



2.1 .3. The p-Ophiuchi filament 

The p-Ophiuchi molecular cloud is a nearby star-forming re- 
gion. The western part of this cloud is bounded by a fila- 
ment structure (Oph-fil hereafter), formin g an edge-on PDR, 



illuminated by th e B2 HD 147889 star dAbergel et al 
Habart et alj|200l . 



1999; 



2.2. IRS observations and data reduction 

The four PDRs were observed using IRS onboard Spitzer, in 
the low resolution (A = 60 - 127) mode. The data were ob- 
tained in the spectral mapping mode with a spatial resolution 
from 3.6 to 10.2 arcseconds), where the spacecraft moves the 
spectrograph's slit before each integration, in order to cover a 
given area of the extended PDR. Each ful l spectral cube wa s 
assembled using the CUBISM software dSmith et al.l 120061) . 
The original Basic Calibrated Data (BCD) files are from the 
S12 pipeline of the Spitzer Science Center. Calibration was 
achieved using the FLUXCON tables provided by the Spitzer 
Science Center together with the retrieved d ata. Additional cor- 
rections for extended sources were applied (ISmith et alj |2004) 
and the remaining bad pixels were removed by hand. In the 
low spectral resolution mode, IRS is divided in two modules: 
the Short wavelengths Low resolution (SL) module, and Long 
wavelengths Low resolution (LL) module. Each module is sub- 
divided into two orders: SL2 (5 to 8.7 pm) and SL1 (7.4 to 
14.5 pm) for SL and LL2 (14 to 21.3 pm) and LL1 (19.5 to 
38 pm) for LL. For Ced 201, both SL and LL modules were 
used, providing the data from 5 to 38 pm. For NGC 7023-NW 
we assembled the SL data from the PID 28 Spitzer program^ 
with the LL2 order data obtained as "bonus" in the SPECPDR 
program, while the LL1 order of the IRS was observing NGC 
7023-E. For NGC 7023-E and Oph-fil, the SPECPDR observa- 
tions were only performed with the LL module, to complement 
the available and good-quality ISOCAM spectral cubes. The 
full spectral cubes, ranging from 5 to 38 pm, were constructed 
assembling the cubes from each module, and reprojecting the 
data on a grid with the lowest resolution (i.e. the LL grid for 
Ced 201 and NGC 7023-NW, and the ISOCAM grid for NGC 
7023-E and Oph-fil). The spectral data from 35 to 38 pm is dis- 
regarded because of its poorer quality. The resulting set of data 
consists of four spectral cubes with wavelengths ranging from 
5 to 35 pm. Each spectral cube is thus a 3-dimensional matrix 
C{p x , p y , A), which defines the variations of the recorded data 
with respect to the wavelength A, for each considered position 
with coordinates (p x , p v ) in the cube. 



2.3. Zodiacal light 

The contribution from zodiacal light to the observations can be- 
come significant when looking at relatively faint objects. This 
is the case of Ced 201. Fortunately, this nebula is situated in 
a region of the sky which is slightly contaminated by zodiacal 
light. We however achieved a background correction for this 
object. This consists in subtracting the spectrum from a nearby 
region of dark sky (including background emission and zodia- 
cal light) to the observations. The other objects show stronger 
emission and are thus less subject to relative contamination by 
zodiacal light (less than 10 % at 20 pm). For these objects no 
correction was applied. 

2.4. MIPS-SED and IRAS data 

We have combined the Infrared-Red Astronomical Satellite 
(IRAS) photometry and Multi-Band Imaging Photometer 
(MIPS, onboard Spitzer) data with the IRS data to build the 
spectral energy distribution of Ced 201, between 5 and 100 pm. 
The IRS spectrum of Ced 201 was obtained by averaging the 
full IRS cube (see Sect. l2.21 >. The IRAS data was retrieved from 
the archive. Because the IRAS beam is larger than the whole re- 
gion observed in the spectral mapping mode with IRS, it is in 
principle not possible to compare the fluxes observed with the 
two instruments. Thus we have made the assumption that the 
emission of Ced 201 is dominated by the region close to the 
illuminating star and covered by IRS observations. Therefore 
the flux calibration of the IRAS points at 12, 25, 60 and 100 
pm points was achieved by scaling the 12 and 25 pm points 
to the IRS-LL data. The MIPS-SED cubes were built stacking 
the BCD files from the S12 pipeline. The MIPS-SED intensity 
spectrum was done by averaging the data over the common re- 
gion with IRS. 

3. Blind Signal Separation Methods 

3.1. Problem overview 

Blind Signal Separation is commonly used to restore a set 
of unknown "source" signals from a set of observed signals 
which are mixtur es of these source signa ls, with unknown mix- 
ture parameters dHvvarinen et alj uOOll) . It has e.g. been ap- 
plied in acoustics for unmixing recordings, or in the biomedical 
field for s eparating mixed el ectromagnetic signals produced by 
the brain (iSajda et alJl2004l) . BSS is most often achieved us- 
ing ICA methods . In particular, an algorithm called FastICA 
( Hy varinen 1999) has proven to be efficient for recovering 
source signals. An alternative class o f methods for achie ving 
BSS is NMF, which was introduced in lLee&Seungl (fl999) and 
then extended by a few authors. 

In the field of Astrophysics, ICA has been sucessfully 
used for sp ectra discriminatio n in infrared spectro-imagery of 
Mars ices dForni et al.l 120051). elimina tion of artifacts in as- 
tronomical images dFunaro et al.l 120031) or extraction of cos- 
mic microwave bac kground signal in Planck simulated data 
(Maino et al. 2002). To our knowledge, NMF has not yet 



2 Spectral Line Diagnostics of Shocks and Photon-Dominated 
Regions, PI Giovanni Fazio. 



been applied to astrophysical problems. However, it has been 
used to separate spectra in other application fields, e.g. for 
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magn etic resonance ch emical shift imaging of the human 
brain (Saida et al. 20041) or for analysing wheat grain spectra 
dGobinet et al]|2005l) 

The simplest version of the BSS problem concerns so- 
called "linear instantaneous" mixtures. It is modelled as fol- 
lows: 



X = AS 



(1) 



where X is an m x n matrix containing n samples of m observed 
signals, A is an m x r mixing matrix and S is an r x n matrix 
containing n samples of r source signals. The observed signal 
samples are considered to be linear combinations of the source 
signal samples (with the same sample index). It is assumed that 
r < m in most investigations, including this paper. The objec- 
tive of BSS algorithms is then to recover the source matrix S 
and/or the mixing matrix A from X. 

The correspondence between the generic BSS data model 
([TJ and the 3 -dimensional spectral cube C(p x ,p y ,A) to be anal- 
ysed in the present paper may be defined as follows. Each ob- 
served signal, consisting of all samples available for this signal, 
corresponds to a row of X in Eq. (Q~|). In this paper, the sample 
index is associated to the wavelength A, and each observed sig- 
nal consists of the overall spectrum recorded for a cube pixel 
(Px>Py)- Moreover, each observed spectrum is a linear combi- 
nation of "source spectra", which are respectively associated to 
each of the (unknown) types of dust particles that contribute to 
the recorded spectral cube. Therefore, the recorded spectra may 
here be expressed according to (03, with unknown combination 
coefficients in A, unknown source spectra in S and an unknown 
number r of source spectra. Sections [3. 2 1 and 1331 introduce the 
two BSS methods we have tested on these data, i.e. FastICA 
and an NMF algorithm. 

3.2. FastICA 

FastICA is a statistical BSS method intended for stationary, 
non-Gau ssian and mutua lly statistically independent random 
signals dHvvar inen Il999h . It is expressed for zero-mean sig- 
nals hereafter. In practice, such signals are obtained by first 
subtracting the sample mean of each observed signal to all sig- 
nal samples in the corresponding row of X in £[). For the sake 
of simplicity, the notation X refers to these zero-mean signals 
hereafter. 

The next step of FastICA consists in applying SVD to the 
covariance matrix of the observed data, i.e. in deriving a matrix 
Z of transformed signals defined as 



Z = MX 



(2) 



where M is selected so that: i) all signals in Z are mutually un- 
corrected, ii) each of these signals has unit power and iii) the 
number of signals in Z is r. The value of f is selected as follows. 
When applied to the m signals in X, SVD intrinsically yields 
m output components. Keeping all these components therefore 
corresponds to selecting r = m. Instead, if r < m, one may 
choose to only keep the r output components which have the 
highest pow ers, with r selected so that r < r < m (see details 
on p. 129 of iHvvarinen et al.l (1200 lb ). This reduces the dimen- 
sionality of the processed data and allows one to combine the 



following two features: i) using f > r still makes it possible to 
recover all source signals from Z and ii) selecting r < m de- 
creases in Z the influence of noise components which exist in 
real recordings X but were not taken into account in the above 
data model (|2). 

The basic version of FastICA then extracts a first source 
signal from the matrix Z. The criterion used to this end consists 
in maximising the non-Gaussianity of an output signal defined 
as a linear combination of the signals in Z. Therefore, denoting 
z a column of Z corresponding to a given sample index, the 
corresponding sample of the output signal reads 



y = d T z, 



(3) 



where the column vector d is constrained to have unit norm. 
Several versions of the FastICA method have been defined, 
depending on which parameter is used to measure the non- 
Gaussianity of y. The most standard parameter is the absolute 
value of the non-normalised kurtosis, defined for a zero-mean 
signal y as 



Kurt(y) 



E[y 4 ] - 3(E[y 2 ]) 2 



(4) 



where E[.] stands for expectation. Various algorithms may then 
be used for adapting d so as to maximise that absolute kurto- 
sis parameter. Before the Fa stICA algorithm was introduced, 
Delf osse & Loubaton (1995) optimised it by using a standard 
gradient ascent procedure. FastICA is an alternative, fix ed- 
point, optimisation algorithm described in lHyvarinenldl999l) .It 
has been shown to yield much faster and more reliable conver- 
gence than gradient procedures. Moreover, it does not require 
one to select any tunable parameter (such as the adaptation gain 
of gradient algorithms). Once a first source signal has thus been 
extracted as an output signal y, FastICA removes its contribu- 
tions from all observed signals contained by X. This yields a 
matrix X' which only contain r — 1 sources. The same proce- 
dure as above is then applied to X' in order to extract another 
source. This "deflation" procedure is repeated until all sources 
are extracted from the observations. 



3.3. Non-Negative Matrix Factorisation 

Unlike ICA, NMF is based on the assumed non-negativeness 
of the source signals and mixing coefficients without requiring 
independence of the source signals. It aims at recovering the r 
source signals by approximating the supposedly non-negative 
matrix X with the following factorisation: 



X&WH 



(5) 



where W is a m x r non-negative weight matrix and H is a r x n 
non-negative matrix of approximated "source" signals. The ap- 
proximation quality in eq. (0 can be optimised by adapting the 
non-negative matrices W and H so as to minimise the squared 
Euclidian distance \\X — Wi/|| 2 or the divergence D(X\WH) 
dLee & Se ung 2001), defined as 



\\X-WHf = Yj(X ij -(WH) i jf 



(6) 
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and 



D(X\WH) = £(X y log -i- - Xtj + {WH)ij) 



{WH)i 



(7) 



The adjustement of W and H can be achieved using classi- 
cal gradient descent but can be problematic (see Lee & Seung 
2001). Therefore, the authors have proposed a new algorithm 
showing that the Euclidian distance is non increasing under the 
iterative update rule 



(W T X) m m m (XH') U 



-,W ia «- W,v 



„„ (WUH 1 ),,, 
and that the divergence is also non increasing under the rule 



(8) 



lZiW ia X ifi /(WH) ifi 



£ M H all X ilA /(WH) it 
v w 



(9) 



Thus, the following iterative algorithm can be derived from 
this result in order to minimise either the euclidian distance or 
divergence: 

1. fix r 

2. initialise randomly matrices W and H 

3. update these matrices with the update rule ([8j or (O 

4. if convergence is reached, then end. Otherwise go back to 
Step 3. 

When convergence is reached, H provides an approximation of 
"source" signals. 

4. Application of the BSS methods 

4.1. Suitability of BSS methods for analysis of 
Spitzer-/F?S cubes 

In order to apply the above BSS methods to the IRS data 
cubes, it is necessary to make sure that the "linear instanta- 
neous" mixture condition is fulfilled (see Sect. I3.ll ). Here we 
consider that each observed spectrum is a linear combination 
of "source spectra", which are due to the emission of differ- 
ent populations of dust particles. The main effect that can dis- 
turb the linearity of the mo del is radiative transfer as shown by 
Nuzillard & Biiaouil feood) . secause of the non-linearity of the 
equations. In our case however, this effect is not expected to be 
a major concern. The mid-IR emission of dust particles occurs 
in the external layers of the clo ud i.e. at optical d epth lower 
or comparable with H2 emission dJoblin et al ] |2007h . Using the 
Meudon PDR code (Le Bourlot et al. 1993, Le Petit et al. 2006) 
we found that H2 emission occurs at extinctions of Ay < 0.2 
in NGC 7023-E and Oph-fil. For Ced 201 which is a partic- 
ular case, as the star has penetrated the cloud, we use the vi- 
sual extinctio n Ay = 0.21 of the illuminating star BD+69° 123 1 



derived bv lWitt et al.l (11987). From this we conclude that very 
small particles emit in the infrared at a maximum visual extinc- 
tion of Ay ~ 0.2 in the studied PDRs. We can then roughly esti- 
mate T9.7, the optical depth of the 9.7 fim silicate band which is 
the main absorption in the mid IR range. Typical ly, the T9.7 /Ay- 
ratio is between ~ 10 and ~ 20 (Mathis 1990l) . which yieds 



T9.7 < 2 x 10 ~ 2 . This value is very small, indicating that the 
effect of radiative transfer on the mid-IR spectra is very weak. 
Thus, the source spectra are not significantly affected by this 
reabsorption, and we conclude that the "linear instantaneous" 
mixture model is fulfilled. This is validated a posteriori as we 
show that the extraction has the same efficiency around 10 fim, 
where the silicates absorb, as at the other wavelengths (see 
FigHO). 

4.2. Preliminary tests 

To test the BSS methods presented in Sect. [3] we generated 
two synthetic spectra: one representing a PAH-type emission 
spectrum, consisting of Lorentzian profiles at 6.2 pm, 7.7 pm, 
8.6 pm and 11.3 /mi, and one representing the emission of 
VSGs (RJB), consisting of a broad emission band at 7.8 pm, 
the 11.3 fim band and a second-order polynomial continuum. 
These spectra were mixed with a random 100 x 2 mixing ma- 
trix A, thus providing 100 artificial observed spectra, which are 
linear combinations of the initial spectra. These spectra were 
then analysed using FastICA and NMF With both methods, 
we managed to recover the original spectra efficiently. In this 
case we have artificially mixed the sources, and therefore we 
can compare the extracted spectra with these sources to mea- 
sure the efficiency of the methods. With both methods, the cor- 
relation between the original and unmixed signals is excellent 
(correlation coefficients higher than 0.995). However a realistic 
test must include the effect of noise. The noise in the IRS in- 
strument is quite complex, and strongly depends on wavelength 
because of the different modules that are used. In order to sim- 
ply quantify the effect of noise on the BSS we applied a white 
spatially homogeneous noise. Though this is not representative 
of the real IRS noise wich would require a detailed study to 
be estimated, this provides good knowledge on the response of 
the algorithms to noise. With both methods the results are still 
good for Signal to Noise Ratios (SNRs) over 2 (i.e. 3 dB). As 
a comparison, the estimated IRS instrumental noise using the 
SpecPET tooJl developed by the Spitzer Science Center, is be- 
tween 20 and 100 (13-20 dB) depending on the wavelength for 
Ced 201. We have also estimated the noise using the spectrum 
of Ced 201 presented on Fig.[9](position 2) and found a SNR of 
about 30 (14.8 dB) in the LL region which is consistent with the 
SpecPET estimation. Below the limit SNR equal to 2, the effi- 
ciency of FastICA drops dramatically while NMF is still able 
to significantly recover the original signals (correlation coeffi- 
cient greater than 0.85). 



4.3. Application to Spitzer data 

We have applied FastICA and NMF to the Spitzer data. The 
spectra from a given cube of observations C(p x ,p y ,A) are 
placed in the rows of the X matrix defined in Eq. ([1) or (0. 
To maintain the efficiency of the methods the noisiest spectra 
were removed. Thus, X contains the m IRS spectra of a con- 
sidered PDR, over n points in wavelength. For NMF, Euclidian 
distance or divergence were minimised using Lee & Seung's al- 



http://ssc.spitzer.caltech.edu/tools/specpet/ 
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gorithm ( Lee & Seung|2001 ) implemented with Matlab, to find 
W and H. The best results were found using the divergence cri- 
terion after about 1000 iterations (which takes less than one 
minute with a a 3.2 GHz processor). The value of r (number 
of "source" spectra) is not imposed by the NMF method. The 
strategy for extracting the sources was the following: 

- Apply the algorithm with a minimum number of assumed 
sources f — 2, to a given dataset, providing 2 sources. 

— > If the found solutions are physically coherent and lin- 
early independent we consider that at least f — 2 sources 
can be extracted. 

— > Else, we consider that the algorithm is not suited for 
analysis (never occurred in our case). 

- try the algorithm on the same dataset but with r — 3 
sources. 

— » If the found solutions are physically coherent and lin- 
early independent we consider that at least r — 3 sources 
can be extracted. 

— > Else, we consider that only two sources can be extracted, 
extraction was over at r — 2 and thus r -2. 

- same as previous step but with f — 4 sources. 

— > If the found solutions are physically coherent and lin- 
early independent we consider that at least r — 4 sources 
can be extracted. 

— » Else we consider that only three sources can be ex- 
tracted, extraction was over at r — 3 and thus r = 3. 

Physically uncoherent spectra exhibit sparse peaks (spikes) 
which cannot be PDR gas lines. We found r — 3 for NGC 7023- 
NW and r — 2 for the other PDRs, implying that we could 
respectively extract 3 and 2 spectra from these data cubes. 
Similarly, with FastICA, all sources were extracted with r — 2 
for Ced 201, Oph-fil and NGC 7023 E, and r = 3 for NGC 
7023-NW. Thus we were able to extract the "source" signals 
from the original cubes using both methods. The choice of the 
method was therefore based on its suitability, considering the 
information we have on the sources. In the following section, 
we argue on why NMF is more consistent in our case. 

4.4. Choice of the BSS method 

NMF and FastICA have different constraints. The strength of 
FastICA is that the algorithm is guaranteed to converge to- 
wards a solution which yields separated spectra. However, this 
property only holds if the source signals are statistically in- 
dependent, due to the non-gaussianity criterion optimised by 
FastICA. This criterion is not ideal in our application, where 
the "source" spectra are likely to be correlated due to chemical 
similarities between the different emitting populations (RJB). 
Moreover, the presence of noise can degrade the efficiency of 
the separation with FastICA (see Sect. 14.21 ). On the contrary, the 
non-negativity of the "sources" and mixing coefficients, used 
in NMF, are consistent with the non-negativity of emission 
spectra and of their linear combination. Thus NMF is likely 
to be more appropriate for the analysis of our data cubes. The 
only identified drawback of NMF is that its convergence point 
may depend on the random initialisation of the matrices (see 




20 
k {\im) 



Fig. 1. BSS extracted spectra in Ced 201. 

Sect. 13.31 1. To test this effect, we have ran NMF on the data of 
NGC 7023-E using 200 different initialisation conditions. We 
found that the outputs only slightly vary with the initialisation 
conditions. Since NMF is most likely to have the greatest effi- 
ciency for our application, we now detail the results obtained 
with this algorithm when applied to our data cubes. 

5. Results 

In this section we present the extracted spectra for Ced 201, 
Oph-fil, NGC 7023-E and 7023-NW (Figs. [BID and their dis- 
tribution maps. These spectra are normalised to 1. This has to 
be done because they have a different intensity depending on 
the considered pixel. The molecular hydrogen lines have been 
substracted in order to provide the spectra of dust only and the 
hydro gen/dust interactio n will be the subject of a subsequent 
paper Ooblin et alj|2007b . 



5.1. Extracted spectra 

From each of the cubes of Ced 201, Oph-fil and NGC 7023- 
E, we were able to extract two spectra. These independently 
found pairs of spectra have the same characteristics from one 
object to another. Moreover, for each object, the two spectra 
exhibit different features. The Signal 1 spectra, presented in the 
upper part of Figs.HE], are characterised by a rising continuum 
combined with broad emission bands located at 6.2, 7.8, and 
11.4 /mi (see Table [B- On the contrary, the Signal 2 spectra, 
presented on the lower parts of Figs. [T]|3] are dominated by the 
AIB emission at 6.2, 7.6, 8.6, 11.3, and 12.7 ym. In NGC 7023- 
NW we could extract three signals from the original data. These 
spectra are presented in Fig. [4] Signal 1 is very similar to Signal 
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Fig. 2. BSS extracted spectra in Oph-fil. 





20 
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Fig. 4. BSS extracted spectra in NGC 7023-NW. In this case, 
the 7.7 pm band of signal 1 was artificially normalised to 0.4 to 
be comparable with the spectra extracted in the other objects. 
The dotted line spectra are those extracted from the ISOCAM 
observations of the over-all NGC 7023 PDR by RJB. 

5.2. Spatial distributions of the extracted spectra 

The observed data in matrix X. may be expressed as the prod- 
uct of two matrices, W and H eq. (5). Therefore, each observed 
spectrum from X is defined as a linear combination of the ex- 
tracted spectra in H. In our case, each IRS spectrum at position 
(p x , p y ) in a given cube can be written as follows: 



20 25 
A (urn) 



Obs(p x ,p y ,A) = 2_ i W(Px>Py)nSn{X) 



(10) 



Fig. 3. BSS extracted spectra in NGC 7023-E. 



1 found in the other PDRs studied above: it exhibits a broad 7.8 
pm band and a steep continuum. Signal 2 and Signal 3 spectra 
are dominated by AIBs but with different relative intensities. 
Signal 2 and Signal 3 are respectively dominated by the 11.3 
and 7.6 pm bands. 



where Obs(p x , p y , A) is the observed spectrum from the pixel 
indexed by (p x ,p y ), S n (A) is the n' h extracted spectrum in H 
and w(p x ,p y )„ are the weights relative to spectrum n. In our 
case n is restricted to a maximum of 2 extracted spectra, ex- 
cept in the case of NGC 7023-NW where n = 3. The NMF 
algorithm extracts the sources up to arbirary scaling factors, 
i.e. it provides y„ = rj p S p , where j] p is an unknown scale fac- 
tor and S p is the p' h "source". Each output can be compared 
to the observations using the correlation between observations 
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and NMF output. We centered the observations (and therefore 
the outputs) and considered that the extracted spectra are not 
correlated. Thus the correlation parameter c p we consider here 
can be written: 

c p = E[Obs(p x , p y , A)y p ] = i] P w(Px, p y ) P E[S P (A) 2 ] (11) 

where E[.] stands for expectation. We can then calculate 
y(p x ,p y ) = Cp/Cp> the ratio of the correlations of Signal p and 
Signal p' to the observations: 

E[Obs(p x ,p y ,A)y p ] 

7(Px,Py) = „ r „, , -^7 T 

E[Obs(p x ,p y ,A)y p >] 

w{p x ,p y ) p E[Sp(A) 2 ]j] p 

= — x ; — - — - (12) 

w(Px,p y ) P > E[S p (A) 2 ]rj p > 

K 

where k is a constant for a given object, y traces the ratio 
of the weights w(p x ,p y ) at the position ( Px , Py ), and thus the 
emission ratio between sources p and p', up to the k scale fac- 
tor, which does not depend on position (p x ,p y ). The maps of 
Figs. |5]|7] show the spatial distribution of the y emission ratio 
between Signal 1 and Signal 2, and Fig.|8]shows the c p corre- 
lation map for each signal. 

In Ced 201 the star is located in the central part of the neb- 
ula, and the Signal 2 emission is also concentrated in this re- 
gion (Fig. [5]), whereas Signal 1 dominates in the periphery of 
the nebula (Pig. O. In Oph-fil, the star is West of the filament. 
Again, Signal 2 is dominant in regions closer to star (Fig. |6]) 
whereas Signal 1 emits in deeper regions (Fig. [6}. In NGC 
7023-E, the illuminating star is located West of the filament. 
The distribution maps show that the maximum Signal 2 emis- 
sion is on the edge of the filament (Fig. [7]). Further East in the 
cloud, the emission is highly dominated by Signal 1 (Fig. [7J. 
Finally, in the case of NGC 7023-NW Signal 1 is dominant 
behind the filament, Signal 2 is dominant on the edge of the 
filament facing the star, and Signal 3 in front of the filament, in 
the closest region to the star (Fig. [8j. 

6. Efficiency of the reconstruction of the 
observations using the extracted spectra 

BSS methods are quite new to astrophysics, and for this rea- 
son they should be used carefully. However, several arguments 
confirm that in our particular case, NMF can be used efficiently. 

Contrary to most approaches, there is absolutely no a pri- 
ori information included in the algorithm, and no subjective 
physical constraints on the source spectra, so that the extracted 
spectra are purely based on the information contained in the 
data. Thus, it is quite convincing that we are able to recover a 
spatial structure without any information on this included in the 
algorithm. In this paper, we show that for the studied PDRs, the 
mid-IR emission spectrum can be fitted everywhere using lin- 
ear combinations of two or three simple spectra extracted from 
the considered PDR. Figure [9] shows the reconstruction of the 
observations taken on three points of the Ced 201 PDR. The 
first position is at the periphery of the PDR, the second one is 
closer to the star and the last one is very close to the star. Using 
linear combinations of the extracted signals (see Sect. [3j we 



2.54e+02 




22:13:30 22:13:24 22:13:21 1. 00e + 02 

Rq (2000.0000) 

Fig. 5. Map of the emission ratios Signal 1 /Signal 2 (upper) and 
Signal 2/Signal 1 (lower) in Ced 201. The 5-35 pm integrated 
mid-IR emission is in contours (contours 1, 2, 3 are respectively 
1, 1.5,2 xl0~ 4 W m 2 sr '). The illuminating star is located at 
the center. 



can reproduce the observations accurately, for the three posi- 
tions. The residuals from this fit are shown in Fig.[l0] The ratio 
of the power of the observed signal to the power the of resid- 
uals is about 30 for position 2 (see Fig. [9j>, which is consistent 
with the value of the SNR we have estimated (see Sect. I4.21 . 
The residuals are at noise level, proving that the reconstruction 
from the extracted spectra of observations is efficient. It is re- 
markable that we are able to show that in some regions of Ced 
201, nearly all the mid-IR emission is due to Signal 1. Also 
note that, even though we chose to use NMF for the reasons 
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Fig. 6. Map of the emission ratios Signal 1 /Signal 2 (upper) 
and Signal 2/Signal 1 (lower) in Oph-fil. The 5-35 yum inte- 
grated mid-IR emission of the filament is in contours (contours 
1, 2, 3, 4 are respectively 3.5, 3.7, 3.9, 4.2 xl(T 4 W trT 2 sr _1 ). 
The illuminating star lies South- West of the filament. 



exposed in Sect. 14.41 FastICA which is a completely different 
method provided similar results, confirming their relevance. 



if 



5:7:48 I 
21:2:18 



21:2:5 
Ra (2000) 
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21:1:51 0.00e + 00 



7. Nature of the carriers of the extracted spectra 

7.1. Continuum-dominated spectrum 

In all the extractions, Signal 1 exhibits a very clear association 
of continuum, and wide AIB emission (Tables Q~||2]i. This par- 
ticul arity was attributed t o the emission of carbonaceous VSGs 
bv lCesarskv et al.1 d2000l) in Ced 201, from ISOCAM observa- 
tions covering a wavelength range from 5 to 17 jjm. The authors 
also proposed that these grains should be found everywhere in 
the interstellar medium. RJB have shown that using mathemati- 
cal analysis, it is possible to identify these grains in other PDRs 
than Ced 201, from ISOCAM observations. They also mention 
that in the emission of this population, the position of the "7.7" 
/im band is shifted towards higher wavelengths which is also 
what we observe (Tables[T]|2]l. When compared to spectra 2 and 
3, the main characteristics of the Signal 1 bands are widths typ- 
ically 50 % higher and a clear red-shift of the 7.7 and 11.3 /mi 
bands. The 8.6 yum band is weak or absent and it is not possible 
to extract the 12.7 //m band from the 12-15 yum plateau. The 



Fig. 7. Map of the emission ratios Signal 1 /Signal 2 (upper) 
and Signal 2/Signal 1 (lower) in NGC 7023-E. The 5-35 /mi 
integrated mid-IR emission of the filament is in contours (con- 
tours 1, 2, 3 are respectively 1, 1.2, 1.4 xl0~ 4 W trT 2 sr~'). The 
illuminating star lies West of the filament. 



15-18 yum region is difficult to study in details because of the 
presence of a strong H2 emission at 17 /zm. 

The IRS-LL module enabled us to build spectral cubes up 
to 35 fim (note that the residual fringes observed in the 20-30 
fim region are instrumental effects and are not due to the algo- 
rithm or signatures of the dust). We have shown that a single 
spectrum, thus a single continuum, can reproduce the contin- 
uum emission all over the Ced 201 PDR (see Sect. [6). This 
indicates that the shape of the continuum does not vary signifi- 
cantly with distance to the star, though its intensity varies. From 
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Table 1. Central wavelengths of the emission bands in the extracted spectra. The vl v integrated intensity ratios of the bands 
compared to the 7.7//m band are written in parentheses. The bands were fitted using lorentzian profiles. 



VSGs PAHs (neutral and cation mixture) 



Bands 


Ced 201 


Oph-fil 


7023-E 


Ced 201 


Oph-fil 


7023-E 


6.2 /im 


6.2 


6.21 


6.27 


6.24 


6.23 


6.24 




(0.06) 


(0.27) 


(0.41) 


(0.35) 


(0.38) 


(0.37) 


7.7 yum 


7.87 


7.77 


7.79 


7.66 


7.72 


7.70 




(1.0) 


(1.0) 


(1.0) 


(1.0) 


(1.0) 


(1.0) 


8.6 yum 




8.54 


8.51 


8.60 


8.57 


8.58 






(0.05) 


(0.05) 


(0.13) 


(0.16) 


(0.21) 


1 1.3 yum 


11.41 


11.40 


11.38 


11.26 


11.29 


11.28 




(0.19) 


(0.13) 


(0.13) 


(0.13) 


(0.26) 


(0.30) 


12.7 /jm 








12.59 


12.66 


12.62 










(0.11) 


(0.21) 


(0.25) 


Cont. 


Yes 


Yes 


Yes 


No/weak 


No/weak 


No/weak 



Table 2. Same as Table 1 but for NGC 7023-NW. 



Bands 




This work 






RJB 




VSGs 


PAHs 


PAHs + 


VSGs 


PAHs 


PAHs + 


6.2 pm 


6.24 


6.22 


6.24 


6.29 


6.27 


6.28 




(0.28) 


(0.34) 


(0.35) 


(0.35) 


(0.47) 


(0.40) 


7.7 yum 


7.74 


7.64 


7.61 


7.82 


7.65 


7.63 




(1.0) 


(1.0) 


(1.0) 


(1.0) 


(1.0) 


(1.0) 


8.6 yum 




8.55 


8.57 




8.58 


8.59 






(0.29) 


(0.27) 




(0.31) 


(0.27) 


1 1.3 yum 


11.36 


11.25 


11.18 


11.36 


11.30 


11.18 




(0.10) 


(0.37) 


(0.17) 


(0.28) 


(0.42) 


(0.12) 


12.7 yum 




12.66 


12.67 




12.76 


12.63 






(0.16) 


(0.08) 




(0.17) 


(0.08) 


Cont. 


Yes 


No/weak 


No/weak 


Yes 


No/weak 


No/weak 



this observation, we can conclude that these grains are most 
likely transiently heated. If this was not the case, their tem- 
perature, and therefore the continuum shape, should strongly 
change as a function of the distance from the star, as it hap- 
pens for classical big grains (BGs). However we cannot com- 
pletely rule out the presence of a fraction of grains at thermal 
equilibrium emitting in the 25-35 yum range. On this basis, and 
because of the strong similarity of our Signal 1 spectrum and 
the ISO VSG spectra in the 5 to 16 yum region, either observed 
( Cesarsky et al. 2000h or extracted mathematically (RJB), we 
propose to attribute Signal 1 in the 5-25 /im range to the emis- 
sion of VSGs which appear to be carbo naceous, and proba- 



bly corresponding to the VSGs used in the lDesert et al. (1990) 



model. In the following sectio n we fit the emissio n of the Ced 
201 PDR using the model of Desert et al. ( 1990h in order to 
provide some clues on the contribution of the two populations 
of grains (VSGs and BGs) to the continuum. 



7.2. Modelling the emission of VSGs versus BGs 

The problem of calculating the infrared emission from 
VSGs in the ISM has been addressed by several authors 



(Desert et al. 1990; Sieben morgen et all 119921; iDraine & Li 



2001; Zub ko et al.ll2004 . One question we would like to ad 
dress is whether BGs contribute to the continuum below 35 /im 
in the studied PDRs. Although this point would deserve a full 
study by itself, we provide here first elements. To estimate the 
emission of BGs, we fitted the total IR emission from Ced 201 
gathering the IRS, MIPS, MIP S-SED and IRAS data, using the 
model of Desert et al. ( 1990b . The input parameters we have 
used for the model are listed in Table [3] T he radiation field 



was calculated using the M eudon PDR code (iLe Bourlot et al 



1993c ILe Petit et alj 12006). The stellar spectru m is from the 



Kurucz library with a temperature of 10500 K dKuruczlll991 ) 
and a radius of 2.26Rq. The distance to the cloud, d, was ad- 
justed to obtain the value of Go =300 at the interface as given by 
Young Owl etal.ld2002l) . This gives d=0.0145pc. This PDR has 



the strongest UV field of the three PDRs for which we have the 
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Fig. 8. Map of the correlation c p of each extracted signal with 
observations, Signal 1 in red, Signal 2 in green and Signal 3 in 
blue. The contours are the Infrared Array Camera (IRAC) emis- 
sion at 8 jum showing the filamentary structure of this PDR. The 
illuminating star is situated in the lower left hand corner. In this 
region a slight artifact correlation of Signal 1 with observations 
is seen due to the presence of a continuum emanating from big 
grains at thermal equilibrium. 




Fig. 10. Residuals from the fit presented in Fig.|9](i.e. observa- 
tion minus reconstruction). 



_ 1 50 - 



~ 300 - 




Fig. 9. Reconstruction of the observations on three pixels of 
the Ced 201 cube, using linear combinations of the extracted 
spectra. Solid lines are the observations and circles show the 
reconstruction. Position 1 is at the periphery 30" North of the 
central star, Position 2 is 8" North of the star and Position 3 is 
very close to the star. 



full 5-35 yum data, and should thus host the hottest BGs, likely 



to emit in the mid-IR continuum. The resulting fit (Fig. fTTb 
shows that the emission from VSGs is dominant up to 50 fim. 
However, because this model is not completely adapted to the 
grain populations we consider, and because of the difficulties 
to calibrate the IRAS data, it is not possible to r ule out BGs 



emissi on over 25 fim. Moreover, as mentioned by Zubkoetal 



(2004), the optical properties of BGs and VSGs are not well 
constrained, implying that their size distribution cannot be de- 
rived precisely. Still, the model provides qualitative informa- 
tion on the emitting population, and we can conclude from the 
fit that the emission in the 25 fim IRAS band is dominated by 
stochastically heated VSGs. 

Table 3. Input par ameters for the fit o f the Ced 201 IR emission 
with the model of Desert et al. dl990h . 



Component 


Relative mass abundance 


a" 






PAHs 


3.410" 4 


3.0 


0.4 


1.2 


VSGs 


5.7 10" 3 


2.6 


1.2 


10 


BGs 


0.994 


2.0 


10 


500 



* a is the exponent in the power law of the size distribution: n(a) oc a ° where a is the size in 
nanometers and n(a) gives the numer density of grains with a radius between a and a + da 



7.3. Pure band spectra: PAHs 

Signal 2 is dominated by the AIBs and is similar to the PAH- 
like emission of RJB. Therefore, we propose that Signal 2 is 
due to free PAH emission. The three extracted PAH spectra of 
Figs. 1-3, are very similar. In this case, the 11-14 fim region 
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representative of c ation emission. This was confirmed by 
Flagey et al.l d2006h : 



XfjX m) 



Fig. 11. Spectral energy distribution (5-100 fim) of the Ced 201 
reflection nebula. The continuous line is the emission modeled 



according to the dust model of lDesert et al.l ( 1990) with the pa- 
rameters from Table [3] 



is dominated by the 11.3 and 12.7 fim bands and no plateau 
such as the one found in the VSG spectra is detected. Here also 
the 15-18 fim region is difficult to study because of the high 
contamination by the H2 line at 17 fim line. However, the spec- 
trum extracted in the case of Ced 201 exhibits a lower 11.3 
fim band than in NGC 7023-E and Oph-fil when compared to 
the intensity of the "7.7" fim band (see Table |2). This varia- 



tion o f the 11.3/7.7 ^m r a tio was interprete d by Jobl in et al 
( 1996); ISloan et all dl999h : iBregman & fern! (|2005|) as an ef- 



fect of variation in the relative abundance of neutral and posi- 
tivel y ionised PAHs (res pectively PAHs and PAHs + hereafter), 
and iFlagev et al. (2006) have shown that the effect of size on 
this ratio was weak. Thus, the PAH extracted spectra in Ced 
201, Oph-fil and NGC 7023-E are mixtures of neutral and 
cationic PAHs. In NGC 7023-NW, RJB were able to identify 
the emission spectra of two populations respectively dominated 
by PAHs and PAHs + . In this object we were also able to ex- 
tract two pure band spectra similar to those determined by RJB 
which we also attribute to PAHs and PAHs + (see FigH|.This 
strengthens the validity of both results. Note however that the 
PAH + spectrum we extract here has a higher 11.3 to 7.7 fim 
integrated intensity ratio than the one of RJB (see Tables [TJ|5J. 
This is because the observations we analyse here do not cover 
the regions close to the star. Thus our PAH + population is "less" 
ionised than the one extracted by RJB, which we will consider 
as reference hereafter. On the other hand, the PAH we have 
extracted is close to the one of RJB (see Table [2] and Fig. @J. 
This extracted spectrum is also similar to th e one observed in 
the HI I region of the Horsehead nebula by ICompiegne et al.l 
(2006), where PAHs are found to be purely neutral because of 
the extreme abundance of electrons. From this we can conclude 
the following: 

- the PAH + spectrum extracted by RJB is dominated by PAH 
cations. In the following, we consider this spectrum to be 



- the PAH spectrum can be attributed to a population con- 
taining only neutral PAHs. 

7.4. lonisation of PAHs in Ced 201, Oph-fil and NGC 
7023-E 



As mentioned in Sect. [73] the PAH spectra of Ced 201, Oph- 
fil and NGC 7023-E are mixtures of neutral and cation PAH 
emission. The reason why we were not able to extract the in- 
dividual PAH and PAH + spectra in these objects is that the 
spectral variations are not strong enough in these cases to pro- 
vide sufficient information to the algorithm to disentangle the 
two populations. In order to estimate the normalised fraction of 
positively ionised PAHs n+/(«o + «+) in these PDRs, we used 
the PAH and PAH + spectra extracted from the data of NGC 
7023-NW by RJB. By comparing the vl v integrated intensity 
ratios between the 7.7 and 11.3 fim bands of these extracted 
spectra (Table 13, to the pure band spectra of Ced 201, NGC 
7023-E and Oph-fil (TableQ]), we could estimate n+/(«o + «+)> 
using NGC 7023-NW as a reference. This ratio is equal to 0.97, 
0.53, and 0.40 for Ced 201, Oph-fil and NGC 7023-E respec- 
tively (Table [3J. This result shows that Ced 201 has a much 
higher proportion of ionised PAHs, probably because of its low 
density (see Table H} i mplying a low r ecombination rate with 
electrons. RJB and Ijoblin et al. ( 19961) also show that a large 
proportion of ionised PAHs is found in NGC 7023-NW and 
NGC 1333-SVS3 respectively. This is likely due to the fact 
that these PDRs are close to the illuminating star. NGC 7023-E 
and Oph-fil, are dominated by neutral PAHs, probably because 
these PDRs are far from the star and have a higher density with 
res pect to Ced 201 (see Table [4]). 

Flag ey et al. ( 20061) have defined the R7.7/11.3 parameter as 



the ratio between the AI V integrated intensity of the 7.7 and 
1 1.3 fim band. We have calculated th is parameter f or Ced 201, 
NGC 7023-E and Oph-fil (Table 151 IFlagev et al.l d2006l) have 
shown that the R7.7/11.3 parameter is related to the ionisation 
parameter Go Vr /n e , where Go is the UV field in units of the 
Habing field, T is the gas temperature and n e the electronic den- 
sity. Using this result, we calculated the ionisation parameter 
from the PAH spectra of our PDRs with the R7.7/11.3 ratio from 
Table|5j We also calculated this parameter independently using 
the physical parameters from Table [4] For this c alculation, we 
used a n e jnu ratio of 1.4 10 -4 dSnow & Witilll995l) . The results 
of these two calculations are reported in Table|5] The estimates 
of the ionisation parameter found here are consistent with the 
ones calculated from the physical parameters. This shows that 
our method to quant ify the ionisa t ion state of PAHs is consis- 
tent with the ones of Flag ey et al. I (l2006l) and could be applied 
to any interstellar PAH spectrum. 



8. Chemical evolution of the very small dust 
particles in the observed PDRs 

The distribution maps presented in Sect. I5.2l clearlv show that 
VSGs and PAHs°/PAHs + are not situated in the same regions 
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Object 


Star: ST, T(K) 


n H (cm 4 ) 


T(K) 


Go (Habing) 


Ref. 


Ced 201 


B9.5V, 10 400 


4.10 2 


200 


300 


Youne Owl et al. (2002) 


Oph-fil 


B2V, 22 000 


3.10 4 


310 


100 


Habart et al. (2003) 


NGC 7023-E 


B5e, 15 000 


7.10 3 


340 


87 


Raoacioli et al. (2006) 



Table 5. Estimated parameters in the studied PDRs: the R7.7/11.3 band intensity ratio, the ionisation parameter, the UV field and 
the ionisation fraction. 



Object 


R7.7/11.3 


Go Vr/n e 


Go 


n + l{n Q + rc + ) 


Ced 201 


3.76 


105000',75000 2 


300 


0.97 3 


Oph-fil 


1.89 


1060',420 2 


100 


0.53 3 


NGC 7023-E 


1.63 


1040\1650 2 


87 


0.40 3 



Using the results from the model described in Flagey et al. i 2006) 2 Using the data from TableRI 



'Using the method described in this paper 



of the studied PDRs. Going towards the star in the PDRs, 
VSGs followed by PAHs and eventually PAHs + (case of NGC 
7023-NW) are successively dominant. This evolution of the 
dust populations across the PDRs while the UV flux is be- 
coming greater suggests that these three populations are chem- 
ically linked. The disappearance of VSG emission, while the 
PAH emission increases, suggests that there is transforma- 
tion of VSGs i nto PAHs unde r the a ction of the UV field, 
as proposed by ICesarsky et alJ (2000) for Ced 201 and by 
RJB for pOph-SR3 and NGC 7023-NW PDRs. Here, we re- 
port new evidence of this transformation in three regions: Ced 
201, NGC 7023-E and Oph-fil PDRs. This is consistent with 
a scenario in which VSGs could be PAH clusters as sug- 
gested by RJB. Interestingly, recent laboratory experiments 
have shown that coro nene clusters can be pho to-evaporated into 
free coronene units (Brech ignac et al]l2005l). PAH aggregates 
would likely form in cold dense regions dBoulanger et al.ll 990: 
Bernard et al. 1993 : Rapacioli et alJ 120061) . and their photo- 
chemistry would start when the UV field is strong enough. 
Thus , the observed var i ations of the mid-IR interstellar spec- 
tra (jPeeters et alJ 120021 : IWerner et al.1 12004 ; Bregman & Temi 
2005!) have to be explained as the evolution of a mixture of 
PAHs and VSGs under the effect of UV flux. The processing 
of VSGs can explain the diminution of contin uum emission 
while a pproaching the irradiating star reported by Werner et al] 
d2004l) . iBregman & Temil d2005l) have found that the 7.7 /im 
emission band shifts towards shorter wavelengths as the UV 
field increases, which they attributed to the evolution from an- 
ionic to cationic species. According to our scenario, this shift 
can be explained as the chemical evolution of the mixture, from 
VSGs, with a band between 7.74 and 7.87 /vm, into PAHs with 
a band between 7.63 and 7.70 yum. In that previous work, the 
authors show that the R7.7/11.3 ratio is increasing when mov- 
ing away from the exciting source, tracing the ionisation state. 



However, after a certain distance, the ratio reaches its maxi- 
mum and starts decreasing again. This was explained by the 
authors as a result of a greater proportion of anionic PAHs, but 
is more likely consistent with the presence of VSGs which ex- 
hibit a low R7.7/11.3 band ratio (see Tables ITT21 . 

9. Conclusion 

Using the data from the Infrared Spectrograph onboard Spitzer, 
combined with powerful Blind Signal Separation methods, we 
were able to extract the spectra of two types of very small 
interstellar dust particles: one carrying mainly the AIB fea- 
tures and the other a mid-IR continuum and broad AIBs. These 
two populations are identified as PAHs and VSGs respectively. 
Concerning the pure AIB spectra, we could identify in the case 
of NGC 7023-NW the emission from PAH and PAH+ dom- 
inated populations. Using the spectra of RJB, we were able 
to simply estimate the ionisation fraction of PAHs in various 
PDRs. Cesarsky et al. (2000) proposed that VSGs should be 
found everywhere in the interstellar medium. However their de- 
tection is far from being obvious. Here we have shown that BSS 
methods enable us to extract their spectrum taking advantage of 
both the spatial and spectral information available in the IRS 
spectral cubes. Thanks to the wide range of wavelengths cov- 
ered, it was possible to confirm that these grains are responsible 
for the interstellar mid-IR continuum and could dominate the 
emission up to 50 //m in cool PDRs. The similarities of their 
spectral features with the AIBs show that they are carbona 



ceous. This is consistent with the predictions of iDesert et al 
dl986l) that VSGs are mostly graphitic and not silicates. The 
distribution maps of PAHs and VSGs further support this con- 
clusion. We indeed show that there is a transformation of VSGs 
into PAHs, probably under the effect of the UV flux. PAHs and 
VSGs clearly dominate the emission in the 12 and 25 fim IRAS 
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ban ds respectively. This result is consistent with the work of 
Fuente et al.l (|1992) who proposed that the 25 fim emission in 
reflection nebulae was due to another type of grains than the 
o nes emitting at 12 um. It is also consistent with the results 
of IVerter et al. (2000), who have shown that the cloud-to-cloud 
variations of the mid-IR emission of molecular cirrus can only 
be explained by changes in the relative abundance of PAHs and 
VSGs. The evolution from VSGs to PAH molecules can explain 
the commonly observed mid-IR spectral variations. The re- 
cently available data from the MIPS photometer in the Spectral 
Energy Distribution mode (MIPS-SED) is currently under anal- 
ysis in order to learn more about the properties of VSGs at 
longer wavelengths. The PAH , PAH + and VSGs spectra ex- 
tracted here can be used as a basis to probe the composition 
of the very small particles in various environments, including 
external galaxies. To a certain extend, this could provide infor- 
mation on the local physical conditions. 
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